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We study the generaUzed quantum isotonic oscillator Hamiltonian given by = -d^/dr^ +1(1 + 
l)/r^ + w^r^ + 2g{r^ — a^)/(r^ + , g > 0. Two approaches are explored. A method for finding 
the quasi-polynomial solutions is presented, and explicit expressions for these polynomials are given, 
along with the conditions on the potential parameters. By using the asymptotic iteration method 
we show how the eigenvalues of this Hamiltonian for arbitrary values of the parameters g, w and a 

may be found to high accuracy. 
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I. INTRODUCTION 

' Recently, Carinena et al. [l[ studied a quantum nonlinear oscillator potential whose Schrodinger equation reads 



o 

(N 



dx^ (2x2 + 1)2 



4>n{x) - E„^Pix) (1) 



The interest in this problem came from the fact that it is exactly solvable, in a sense that the exact eigenenergies and 
eigenfunctions can be obtained explicitly. Indeed, Carinena et al. [l[ were able to show that 

^1 ■ I Vn[X) - (2j;i + i)e 



(N 

[E„^-3 + 2n, n = 0,3,4,5,. 
where the polynomials factors P„(x) are related to the Hermite polynomials by means of 



(2) 



P (x) = l ^ if n = , . 

"^^ \ +4nff„_2(x) +4n(n-3)i/„_4(x) if n ^ 3, 4, 5, . . . ^''^ 

In a more recent work. Fellows and Smith 6] showed that the potential V{x) = x^ + 8(2a;2 — l)/(2a;2 + 1)2 -^^rell as, 
for certain values of the parameters w, g and a, the potential V{x) = w'^x'^ + 2g(x'^ — a?)/{x'^ + a2)2 of the Schrodinger 
equation 

r J2 _2 „2 1 

^n{x) = 2En4'{x), (4) 



+2g (,2+,2)2 



are indeed supersymmetric partners of the harmonic oscillator potential. Using the supersymmetric approach, the 
authors were able to construct an infinite set of exact soluble potentials, along with their eigenfunctions and eigen- 
values. Very recently, Sesma [ol], using a Mobius transformation, was able to transform Eq.([5]) into a confluent Heun 
equation [8| and thereby obtain an efficient algorithm to solve the Schrodinger equation @ numerically. 



'Electronic address: nsaad@upe i.ca| 



t Electronic address: rliall@matlista t.concordia.cal 
t Electronic address: liciftci@gazi.edu . tr| 
§ Electronic address: ,yesiltas@gazi.edu.tr| 



2 



The purpose of the present work is to provide a detailed sohition, by means of the quasi-polynomial solutions and the 
application of the asymptotic iteration method [H^^ j for the Schrodinger equation 



l{l + l) 



+ w^r"^ + 2g 



Tp{r) = 2E'iJj{r), 



(5) 



where I is the angular momentum number Z = — 1, 0, 1, . . . . Our results show that the quasi-exact solutions of Sesma 
as well the results of Carifiena et al. ,1] follow as special cases of our general approach. The present article is 
organized as follows. In the next section, some preliminary analysis of the Schrodinger equation ([5]) is presented. A 
general approach for finding polynomial solutions of Eq.Q, for certain values of parameters w and g, is presented, 
and is based on a recent work of Ciftci et al. [3] for solving the second-order linear differential equation 



a2,ix y 



(6) 



More general quasi-exact solutions, including the results of Sesma are discussed in section III. Unrestricted 
solutions of Eq.Q based on the asymptotic iteration method are discussed in Section IV. 

II. GENERALIZED QUANTUM ISOTONIC OSCILLATOR - PRELIMINARY RESULTS 

A simple scaling argument, using r — a^x, allows us to write the equation ([5]) as 



+ {wa'^fx'^ + 2g 



(x2 + 1)2 



■0(a;) = 2Ea^'il;{x) 



(7) 



A further substitution z — 



1 yields a differential equation with two regular singular points aX z — Q,l and one 



irregular singular point of rank 2 at z = oo. The roots /z's of the indicial equation for the regular singular point z — 
reads ^± = |(1 ± VI + 4.g), while the roots of the indicial equation at z = 1 are ^+ = [l + l)/2 and /i_ = —1/2. 
Since the singularity for z — > oo corresponds to that for a; — >■ cxd, it is necessary that the solution for z — )• oo behave 
as ^^{x) ~ eyi]}{—wa?'x'^/2). Consequently, we may assume the general solution of equation ([7|) which vanishes at the 
origin and at infinity takes the form 



^^{x) = x'+'{x' + lYe-—- Ux). 



(8) 



A straightforward calculation shows that fn{x) are the solutions of the second-order homogeneous linear differential 
equation 

2{l + 1) A^lx 



2£;a2-u;a^(2/ + 3 + 4/i) 



2^(2/ + 3 + 2wa2) + 4^(^ - 1) - 2g 4(.g - fi{p - 1)) 



1 



(2)2 + 1)2 



fix) = 0. (9) 



In the next sections, we attempt to give a general solution of this equation. For now, we assume that fi takes the 
value of the indicial root 



which allows us to write Eq.® as 



/"(-)+ ('^^ + ^-2-«^- )/«(-) + 



2Ea'^ -wa^{2l + 3 + 4fi) 



2fi(2l + 3 + 2wa2) + 2fi{p - 1) 



(10) 



fnix) = 0. 
(11) 



We now consider the cases where the following two equations are satisfied 

2^{2l + 3 + 2u;a2) + 2fi{fi - 1) = 0, 

5 = - 1)- 
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The solutions of this system, for g and /x, are given exphcitly by 

g = 0, ( g = 2{\ + l + a'^w)(i + 21 + 2a'^w), 

or < 

/i = 0, [ -2(1 + ; + a^w). 

In the next, we consider each case of these two sets of solutions. 



(12) 



A. Case 1 

The first set of solutions (g, ^) = (0, 0) reduces the differential equation ^ to 

xf'ni^) + [-2wa'^x^ + 2{l + l)]f^{x) + {2Ea'^ - wa^{2l + 3)) x fn{x) = 
which is a special case of the general differential equation 

(03,02^^ + a3,i^^ + ^3,2X + 03^3) y" + {a2.oX^ + 0342; + a2.2) ?/ - {rifiX + rij) y = 0, 



(13) 



(14) 



with 03,0 = 03,1 = 03,3 = a2,i = Ti^i = 0, 03,2 = 1, 02,0 = —2wa^, 02,2 = 2{l + 1), and n^o = —2Ea'^ + wa'^{2l + 3). 
The necessary and sufficient conditions for polynomial solutions of Eq. (jl4p are given by the following theorem 



Theorem 1. The second-order linear differential equation |j^|) has a polynomial solution of degree n if 



Ti,o = n{n - l)a3,o + na2,o, n = 0, 1, 2, . . 
along with the vanishing of (n + 1) x (n + 1) -determinant A„+i given by 

l3o ai 771 
71 Pi a-2 772 

72 P2 013 r]3 



(15) 



where 



ln-2 13,1-2 Oin-l Vn-1 
7n-l /3n-l On 
7n Pn 



Pn = n,i - n{{n - 1)034 + 0-2,1) 
Cin -n{{n - l)a3_2 + 02,2) 
In = Ti,o - {n - 1)(('T- - 2)a3,o + 02,0) 
Tin = -n{n + 1)03,3 



= 



(16) 



and Ti o is fixed for a given n in the determinant A„+i = 0. 

Thus, the necessary condition for the differential equation (jl3p to have polynomial solutions fn{x) = "Yll^^o^i^^ 

2Ena^ ^wa^{2n' + 21 + ?,), n' = 0,1,2,... (17) 
while the sufficient condition, Eg p^ . is 



A. 



n+l 



ai 

71 q;2 

72 as 



7n-2 



7n-l 

7 



a„_i 








if ri = 0,2, 4,... 



n (-l)2J+ia2j+i72j+i =0 if ji = 1, 3, 5, . 
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where /?„ = 0, = —nin + 2Z + 1) and 7„ = 2wa?'{n — n' — \). 



li I = —1, the determinant A„_(_i is identically zero for all n, which is equivalent to the exact solutions of the 
one-dimensional harmonic oscillator problem. 

For I ^ —1, we have for 7i = 0, 2, 4, . . . , A„+i = and we obtain the exact solutions of the Gol'dman and Krivchenkov 
(or Isotonic) Hamiltonian Hq where 



2 4 2 

wax 



■ipnlix) = 2E'^i°a'^'lpnl{x), < X < CO. 



dx'^ x^ 
These exact solutions are given by 

2a'^El=° = wa^{^n + 2Z + 3), n = 0, 1, 2 



(18) 



(19) 



^„;(x) =a;'+ie-""^ ^ /^i^i(-n; / + |; wa^x^), n = 0,1,2,.... 

where the confluent hypergeometric function i_Fi(— n; a; z) defined, in terms of the Pochhammer symbol (or Gamma 
function r(a)) 

r(a + fc)_/l if (fc = 0,aeC\{0}) 



as 



r(a) \ a(a + l)(a + 2) ... (a + A: - 1) if (fc = N, a e C) 
iFi(-n;a;z) 2^ ■ 



k=0 



{a)kkl 



(20) 



The polynomial solutions fn{x) = iFi{—n;l + ^; waPx"^) are easily obtained by using the asymptotic iteration method 
(AIM), which is summarized by means of the following theorem. 

Theorem 2: (H. Ciftci et al.Q, equations (2.13)-(2.14)) Given Aq = Ao(a;) and sq = sq{x) in C°° , the differential 
equation 

f"{x)^Xo{x)f'{x)+so{x)f{x) 

has the general solution 

C2 + C\ I exp I / (Ao(t) + 2a(T)) dr I dt 



f{x) = cxp — / a{t)dt 



if for some n e N+ = {1, 2, . . . } 

An A„_i 

where 



For the differential equation p3| with 



= a{x), or (5n(a;) = A„Sn_i - An_is„ = 0, 

An = -^n-l + + AgAn, 

Sn — Sn— 1 S^Xn- 

Ao(:r) = - (-^-'";+^<'+^» , 

sa{x) = -{2Ea^ - wa'^(2l + 3), 
the first few iterations with 5n — AnSn-i — An-iSn = 0, using (j2ip . implies 

/o(a;) = 1 

fi{x) = 2u)a2a;2 - (2Z + 3) 

/a (a;) = 4u)2a4a;'* - 4u;a2(2/ + b)x^ + (2Z + 3)(2; + 5) 



(21) 



(22) 



(23) 



(24) 



(25) 



which we may easily generalized using the definition of the confluent hypergeometric function, Ea(|20p. as 



fn{x) = iFi(-n;Z + -;wa^a;^) 



up to a constant. 



(26) 
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B. Case 2 

The second set of solutions 

{g, n) ^ (2(1 + I + a^w){3 + 21 + 2a^w), -2(1 + I + a^w)) 
allow us to write the differential equation ^ as 
'2(/ + l) 8{l + l + a'^w)x 



m + 



„2 



- — '— - 2wa^x j /;(x) + {2Ea^ + wa^{6l + 5 + Swa^)) fn{x) = 0. 



(27) 



A further change of variable z — x"^ + 1 allows us to write the differential equation <\27\ as 

Az{z - l)f"{z) - {Aa^wz"^ + 2(6/ + 5 + Qwa'^)z - + 1 + wa^)) f'{z) + {2Ea^ + wa^{6l + 5 + 8wa'^))z f{z) = 0, 

(28) 

Again, Eq.((28|) is a special case of the differential equation (fT4|) with 03,0 = 03,3 — ri.i — 0, a3_i — 4, 03,2 = —4, 02,0 — 
—A-wa^,ai,\ = — 2(6Z + 5 + 6wa^), 02,2 = 16(/ + 1 + wa^) and ti^q = ~2Ea^ — wa^(6/ + 5 + 8wa^). Consequently, the 
polynomial solutions /„(x) of (l28l) are subject to the following two conditions: the necessary condition reads 



2£;„a^ = 'ii;a2(4n' - 6/ - 5 - Swa^), n' = 0, 1, 2, . . . 
and the sufficient condition; namely, the vanishing of the tridiagonal determinant Eq (|16p . reads 

71 /?! "2 

72 /32 "3 



(29) 



where 



/3„ 



7n-2 /3n-2 

7n-l "n 
7n /3n 



-2n{2n~6l - 7 - Gwa^) 
An{n - 4/ - 5 - 4a^u;) 



= 



7„ = 4ti;a^(n — n' — 1) 



(30) 



and n' = n is fixed for the given dimension of the determinant A„+i. From the sufficient condition pop we obtain 
the following conditions on the parameters 

A2 = ^ a^w{l + 1 + a^w) = 

A3 = ^ a^wH + 1 + a^w)(l + 21 + 20^10) = 

A4 = ^ a^w{l + 1 + a^w){l + 21 + 2a^w){3{l + 61) + Ua'^w) = 

A5 = =^ a^w{l + 1 + a^w){l + 21 + 2a^w){3{6l - 1)(6/ + 1) + 4(38/ + l)a^w + 44aW) = 

Ae = ^ a^w{l + 1 + a^w){l + 21 + 2a^w){3{2l - 1){61 - 1){61 + 1) + 2{208f - 54/ - 5)a^w + 200/aW) = 



For a physically meaningful solution we must have a^w > 0. This is possible for a very restricted value of the angular 
momentum number /. Since /3o = 0, we may observe that 

73 l^s cti 
74 134 as 



A„+i ={l + l + a'^w){l + 2l + 2a'^w)x 



ln-2 Pn-2 

7n — 1 /^n— 1 

7n i^n 



= (/ + 1 + a'^w){l + 21 + 2a'^w) x Q^^_^{a^w) 



6 



where Ql^_i{a'^'w) are polynomials in the parameter product a^w. 

For physically acceptable solutions, we must have I — —1 and the factor {I + 1 + a^w) yields a^w = 0, which is not 
physically acceptable; so we ignore it. The second factor (1 + 2/ + 2a?w) implies a special value of a^w — 1/2, for all 
n, which we will study shortly in full detail. Meanwhile, the polynomials Q^^{a?w) 



1 if n = 2 

lAa^w - 15 if n = 3 

44a''w2 - U%a^w + 105 if n = 4 
200a'*w;2 - hUa'^w + 315 if n = 5 



(31) 



give new values, not reported before, of a^w that yield quasi-exact solutions of the Schrodinger equation (with one 
eigenstate) 



{x^ - 1) 



where 



(a;2 + 1)2 

V'„(x) = (x2 + l)-2-^-e-«"^'-'/2/„(x) 



i^nix) — wa2(4n + 1 — 8a'^w)^pnix) 



For example, A4 = implies, using (I5T]) . that a'^w — y|, and thus we have for 



-ij'^ix) 



225 n 660 (x^ - 1) 



196 49 (a;2 + 1)2 



I ( ^ 465 



the exact solution 

^^(x) = {x^ + l)"'*'e"a"^'(452:'' + 2253;" + 315x'-' - 49) 
with a plot of the wave function and potential given in Figure 1. 




FIG. 1: Plot of the unnormalized wave function ip3{x) and the potential V3 = fffa;^ + ^ (x'^+iyi 



(32) 



and fn{x) are the solutions of 

4z(z - l)/"(z) - (4a2u;z2 + 2(-l + 6wa2)z - 16wa^) f'{z) + inwa^z f{z) = 0, z = x'^ + 1. (33) 



(34) 



Further, A5 = 0, Eq.dSTj) implies 



2 37 V214 
w = — ± „^ 
22 22 
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and we have for 

-i'Ux) + 



^37 ^ V2U 37 V2U A8 7214 jx' - 1) 

Uo =^ on ) ^^'-ii 11 All =C 11 h^2j_T\2 



22 22 



11 



11 "11 



11 '(x2 + l)2 



, , , ,37^ ^/2T4, .39 47214, , . . 



(35) 



the exact solutions 

, + / \ / 9 , , I- 37 _|_ n/2T4 |-37_l_ v/2T4 2 

(a;2 + 1) 11 44 )^ 

(1575a;^ + (9660 ± 420V214)a;^ + (26250 ± 2100V214)a;'* + (29820 ± 2940V2T4)x2 - (1129 ± 188V2T4)). 
Similar results can be obtained for A„+i = 0, for n > 5. 



C. Exactly solvable quantum isotonic nonlinear oscillator 

As mentioned above, for I = —1 and a^w — 1/2, it clear that A„+i = for all n and the one-dimensional Schrodinger 
equation 



.x2 4(.x2 - 1) 
"d!^ ^ T ^ (x2 + 1)2 



^prl{x) = {2n- -)ipn{x), 0,1,2,. 



(36) 



has the exact solutions 

^|;4x) = {x' + l)-'e-^'/'Mx), (37) 
where fn{x) are the polynomial solutions of the following second-order linear differential equation (z — x'^ + 1) 

4z(z ~ l)/;'(z) - (2z2 + 4z - 8) /;(z) + 2nz /„(z) = 0, (38) 
By using AIM (Theorem 2, Eq. ipT]) ). wc find that the polynomial solutions fn{x) of Ea. ([55| are given explicitly as 



(39) 



a set of polynomial solutions that can be generated using 

foix) = 1, Ux) = -3x(2n + l)iFi(-n; | ^{x - 1)) + 6((n + l)x - l)iFi(-n + 1; | ^{x ~ 1)), (40) 

up to a constant factor, where, again, iFi refers to the confluent hypergeometric function defined by (PO)) . Note that 
the polynomials fn{x) in equation (|40| can be expressed in terms of the associated Laguerre polynomials [10| as 



r foix) 


= 1 








hix) 


= a;2 


- 2 






, hix) 


= a;3 


- 6a;2 -1- 
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hix) 




- 16x^ - 


h 52a;2 _ 


52 


hix) 


— 


- 30x4 _ 


h 250x^ 


- 580a;2 + 464 



/o(a;) = 1, /„(x) 



3(-l)"0Fr(n) 
2r(n-f f) 



iil + n)ix - 1)^ +n)a 



x-1 



{x-l)iil + n)x~l)L^ 



s. / X — 1 



(41) 



III. QUASI-POLYNOMIAL SOLUTIONS OF THE GENERALIZED QUANTUM ISOTONIC 

OSCILLATOR 

In this section we study the quasi-polynomial solutions of the differential equation ([SJ. We note first, using the change 
of variable z = x'^, Eq.® can be written as 



2Ea'^ - wa^ (2^ + 3 + 4/i) /i(2; -I- 3 -f 2u;a2) g {z - I) nin - I) 



4z 



2ziz + l) 



2z(z-Hl)2 (2-hl)2 



hiz) = 



(42) 
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By means of the Mobius transformation z = — t) that maps the singular points {— 1,0, cx)} into {0,1, oo}, we 
obtain 



m + 



2/ + 3 2(/^ - 1) 



2t{i-t) 1-t {i-ty 

where we assume 

2Ea'^ -{21 + 3 + 4fi)wa'^ = 
The differential equation can be written as 



(1(21 + 3 + 21110^) 5(2<-l) /i(M-l) 
2i(l-i)2 2t{l-ty ^ (1-^)2 



/„(t)=0, (43) 



(44) 



(t'-2f+t)f:{t) + 



7, 



-2{fj. - l)r + {2fj. -wa^ -I - -)t +{1 + -) 



iiiiii - 1) - g)t + I + ^(/ + ^ + wo") 



fn{t) - 



(45) 



which we may now compare with equation (jl4l) in Theorem 1 with 03.0 = 1,03.1 = —2,03^2 — 1,13.3 — 0, 02,0 — 
-2(/Li - 1), 02,1 = {2fi -wa^-l- 7/2), 02,2 = {I + 3/2), n^o = -(m(m - 1) - .9), n,i = -f - m(^ + | + wa^)- We, thus, 
conclude that the quasi-polynomial solutions fn{t) of Ea. (|45|) are subject to the following conditions: 

g={li-k){(i-k-l), fc = 0,l,2,... (46) 

along with the vanishing of the tridiagonal determinant A„+i = 

71 /3l "2 

72 P2 Olz 

ln-1 Pn-1 OLn 
7n Pn 

where 

Pn = --\{g + {(i-n){3 + 21 + An + 20:^111)), 

an ^ -n{n + l + \), (47) 

7n = g - (m- + i)(Ai-- '^), 

Here, again, g = (/x — fc)(/i — fc — 1) is fixed for given k = n, the fixed size of the determinant A„+i. 



A. Particular Case: n = 

For k {fixed) = n = 0, the differential equation (|45|) has the exact solution /o(i) = 1 if .9 and /i satisfies, simultaneously, 
the following system of equations 

g + n{3 + 2l + 2a'^w) =0, g = (i{(i-l). 
Solving this system of equations for g and fi, we obtain the following values of 

g = 2{l + l + a'^w){3 + 21 + 20^-0)), and (i ^ -2{l + 1 + wa^), 
and the ground-state energy, in this case, is given by Eg.ip^. namely, 

1 



(48) 



Ea^ = --a^w{5 + 61 + Sa^w) 



(49) 



which in complete agreement with the results of Section II. B. 
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B. Particular Case: n = 1 

For k (fixed) = n = 1, the determinant A2 = of (|47)) yields 

g2 + + lo^i + 21(2^1 + 1) + 2a2w(2^ - 1)) + fi{^ - 1)(15 + + 8l{2 + a^w) + 4:a^w{5 + a^w)) = 0, 

(50) 

g-(M-l)(M-2) = 
where the energy is given by use of Ea. (|44l) . for the computed values of /j, and g, by 

E = {l+'^ + 2fi)w. (51) 

Further, Eg. dSOl) yields the solutions for / as functions of /i and a^w 



_ 2 - (5 + 4a^w)/i ~ 2^^ ± ^4 - 4(3 + Sa^w)/! + 9^^ 

^- >-l, (52) 

where the energy states are now given by ([5T|) along with I given by Eg. ([52]). We may also note that for 

a'^w=^{k + l), fc = 0,l,2,... (53) 



and 



1 



a^Ei = (fc + 1) ( -2 + (2fc + l)fj. - ± ^4-4(4A: + 7)Ai + 9^2 _ (54) 

8^ V / 

Further, for g = (fi — 1)(m ~ 2), we obtain the un-normalized wave function (see Eq.®) 

, N , Ox 1 ,n^2„2 /o , 1 + 2/ + /i + 2a'^u) 

^.,(.) = .'-(1 + .^r^^e- l\X + ^^^^^^^^^) (55) 

Thus, we may summarize these results as follows. The exact solutions of the Schrodinger equation ([71) are given by 
Eas. ([5^ and (|55p only if g and fj, are the solutions of the system given by Eg . ((50)) . In Tables HI and HIl we report few 
quasi-exact solutions that can be obtained using this approach. 

C. Particular Case n = 2 

For k (fixed) = n — 2, the determinant A3 — along with the necessary condition (j47p yields 

g^ + 3^2(7/^ - 1 + 2/(1 + /i) + 2a^w(^i - 1)) - .g[18 + 56/ + 8/^ + 18(7 + 2l)^i - 3(5 + 2/)(7 + 2l)^i'^ 

-12a'^w(p - 1)((7 + 2l)fi - 4) - 4a'*w2(2 + 3(p - 2)^)] + fi(p - 2)(p - 1)(105 + 142/ + 60/^ + 8/^ + 6a^w(5 + 2/)(7 + 21) 
+12aW(7 + 21) + 8a^w^) = 0, 

[ g-(M-2)(Ai-3)=0 

(56) 

where, again, the energy is given, for the computed values of /i and g using Egs. (|44p and (j56l) . by 

E= (l+'^ + 2fi)w. 

In Table lllli we report the numerical results for some of the exact solutions of /j, and g using Eq. (f56)) and the values 
of (/, wa^) = (-1, i), (/, wa"^) = (-1, 1), (/, wa^) = (-1, |), (/, wa'^) = (-1, 2) (/, wa^) = (0, i), and (/, wa^) = (0, 2), 

respectively. We have also computed the corresponding eigenvalues E^'f = El^f (fJ-Tg)- 
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TABLE L Conditions on the value of the parameters g and fj, for the quasi-polynomial solutions in the case of A2 = with 
different values of wa^ and I. 



n I wa? Conditions 



1 -1 h 



(^-3 - 15^-1/=^ - yl^/^j , A = 3(36 - %/96T) 
g = 1^-2/3(15 + 6^1/3 + yl2/3)(i5 + 9^1/3 _^ ^2/3) 



M = 



1^-2/3(19 + 8A^/3 ^2/3)(^9 _|_ ^^^1/3 ^ ^2/3) 



El-. 



,(^+2^1/3 + ^^-1/3) 



= i (-7-25vl-^/3„^i/3 



A = 199 - 18^ 
g = iyl'^''3(25 + lOA^/^ ^ ^2/3)^25 + IZA^^'^ + ^^''^) 



1,-1 



,(f + 1^1/3 + 



l/3\ 



= i (^-9 - 33^-'''3 _ ^1/3^ ^ A = 3(72 - yugT) 
g = 1^-2/3(33 + 12^1/3 + yl2/3)(33 + 15^1/3 _^ ^2/3) 



I 



^ = 
= 2 



^^0 = I"' 



M = -i(7+^/T7) 
3 = 29 + 5^17 



£^'o = -l(ll + 2^)t 



M = -i(7-^/T7) 
3 = 29 - 5^17 



g^{-4 + B){~5 + B) 
B = |K (^A^/3 ^ 33^-1/3^ 



-108 + 3i\/2697 



£■1,0 



2B) 10 



/i = -3 - 
g = (5 + i3) (4 + B)) 



B = 3? 



ll(l+iy3)A-l/3 , (1-»x/3)a1/^ 



Si,o = -(| + 2B)ix; 



A : 



-108 + 3i\/2697 



Ai = -3- B, 

<;=(5 + S) (4 + B)) 

^ ^ Sfj ^n^i-n/3)A;^ _l_ (i+»y3)Ai/^ ^ 



A : 



-108 + 3i^2697 



IV. NUMERICAL COMPUTATION BY USE OF THE ASYMPTOTIC ITERATION METHOD 



For the potential paranieters w, a? and g, not necessarily obeying the conditions for quasi-polynomial solutions dis- 
cussed in the previous sections, the asymptotic iteration method can be employed to compute the eigenvalues of 
Schrodinger equation ([7]) for arbitrary values w, c? and g. The functions Aq and sq, using Eas. (j43l) and (|44l) . are given 

by 



Ao(i) = 
so(t) = - 



2i-|-3 
2t(l-t) 



(1-t) 
2t(l-t)^ 



g_(2t-l) 
2 t(\-tY 



(57) 
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TABLE II: Conditions on the value of the parameters g and /i for the quasi-polynomial solutions in the case of A2 = with 
different values of wo? and I. 



wa Conditions 



+ B 



{35-12B)w 



+ 9iV863 



11 

3 



p = i(17 + 3B)(14 + 3i3)) 

B=l?R (43(1 + + (1 - iV3)A^/'' 



Elo = -U^5 + 12B)w 



A = -98 + 9iV863 



i(17 + 3B)(14 + 3B)) 



B = (43(1 - i^/3)yl-l/=' + (1 + iV3)A^^^) , A = -98 + 9i^/863 



Si\o = -s(35 + 12B)w 



2 <; 3= i(-16 + 3B)(-19 + 3B) 



yl = -55 + 165i\/6 



^1*0 = -1(43- 125) w 



M = 
9 = 



i(16 + 3B)(19 + 3B)) 
= 13? (55(1 + iVS)A~^^'' + (1 - iV3)A^^'' 



Elo = -i(43 + 12B)w 



A = -55 + 165iV6 



B, 



3= i(16 + 3B)(19 + 3B)) 

B = |K (55(1 - iV3)yl-i/'' + (1 + iV3)A^/^) , A = -55 + 165i\/6 



Elo = -i(43 + 12B)«; 



where t S (0, 1). The AIM sequence A„(x) and Sn{x) can be calculated iteratively using the iterative sequences ([23]) . 
The energy eigenvalues of the quantum nonlinear isotonic potential ([7]) are obtained from the roots of the termination 
condition (j22p . According to the asymptotic iteration method, in particular the study of Brodie et al , unless the 
differential equation is exactly solvable, the termination condition (1221) produces for each iteration an expression that 
depends on both t and E (for given values of the parameters wa^, g and I). In such a case, one faces the problem 
of finding the best possible starting value t = to that stabilizes the AIM process 0. Fortunately, since t e (0,1), 
the starting value Iq doesn't represent a serious issue in our eigenvalue calculation using (|57p and the termination 
condition p2|) in contrast to the case of computing the eigenvalues using Ao(a;) and 50(2;) as given by, for example, 
equation (jHl), where x G (0, 00). In Table HVl we report our numerical results for energies of the four lowest states 
of the generalized isotonic oscillator of parameters w and a such that wa'^ = 2 and for different values oi g. In this 
table, we set I = —1 for computing the energies Equ^ and £'20^, while we put / — for computing the energies Eio^ 
and E^a'^, respectively. For most of these values, the starting value of t is to — 0.5 and is shifted towards zero as 
g gets larger in value. For the values of g that admit a quasi-polynomial solution, the number of iteration doesn't 
exceed three. For most of the other values of g, the total number of iteration didn't exceed 65. We found that for 
wa^ = 2 and the values of g reported in Table IIVI the number of iteration is relatively small compared to the case of 
wa^ — 1/2 and a large value of the parameter g. The numerical computations in the present work were done using 
Maple version 13 running on an IBM architecture personal computer in a high-precision environment. In order to 
accelerate our computation we have written our own code for a root-finding algorithm instead of using the default 
procedure Solve of Maple 13. These numerical results are accurate to the number of decimals reported. 
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TABLE III: Exact eigenvalues for different values of I and wa^ in the case A3 = 0. 



n I 




Conditions 




^EZf{n,g) 


2 — 1 


1 
2 


gi = 77.22293097048609 


1 

2,- 


— _fi 051 S7nS7SQQ41 QS 






/X2 = —2.4855365082108594 
32 = 24.605574274703333 


i 

771 2 


1 = —2.2355365082108594 




1 


/Ji = -7.398182984326876 
gi = 97.7240263912181 




1 = -7.148182984326876 






^2 — — o.ooouo ( yui^yuoiy^ 
ff2 = 34.03170302988033 


-C'2,- 


— — ^ 1 n^n^7Qni 4Qfi81 Q4 






53 = 2.1530873564462514 


-C'2,- 


— 1 1 QQS1 n^41 7^747^f; 




3 
2 


pi = -8.469623341124414 
ffi = 120.08263624614156 


a 


1 = -8.219623341124414 






H2 = —4.27750521216504 
g-2 = 45.684576900924284 


Ell 


1 = —4.02750521216504 






Hz = 0.9282653601757613 
g-i = 2.2203497780234294 


E2- 


1 = 1.1782653601757613 




2 


m = -9.525122115065386 
gi = 144.35356188223463 


El- 


1 = -9.275122115065383 






fi2 = —5.226942179911145 
g2 = 59.45563545168999 


ir2 

E2- 


1 = —4.976942179911145 






1^3 = 0.9186508169859244 
gs = 2.250665238619284 


I7i2 


1 = 1.1686508169859244 


2 


1 
2 


m = -8.032243023438463 
gi = 110.67814310476818 




1 = -7.282243023438463 






/X2 = -4.32825470612182 

g2 = 46.37506233167478 


S2,- 


1 = -3.57825470612182 




2 


Hi = —11.307737259773461 
gi = 190.4036082349363 


172 

E2- 


.1 = —10.557737259773461 






M2 = -7.180564905703867 

52 = 93.46333689354533 


El- 


1 = -6.430564905703867 






/i3 = 0.9472009101393033 
gs = 2.1611850134722084 


El- 


1 = 1.6972009101393033 



V. CONCLUSION 

We have provided a detailed solution of the eigenproblem posed by Schrodigcr's equation with a generalized nonlinear 
isotonic oscillator potential. We have presented a method for computing the quasi-polynomial solutions in cases where 
the potential parameters satisfy certain conditions. In other more general cases we have used the asymptotic iteration 
method to find accurate numerical solutions for arbitrary values of the potential parameters w and a. 
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TABLE IV: Energies of the four lowest states of the generahzed isotonic oscillator of parameters w and a given for Z = —1 as 
wa^ = 2 and for different values of the parameter g. The subscript mimbcrs represents the number of iterations used by AIM. 



wa^ g 


Eoa^ 




E2a' 


Esa^ 


2 0.000 01 


0.999 993 709 536(39) 


2.999 997 742 768(26) 


4.999 998 464 613(32) 


6.999 998 987 906(23) 


0.1 


0.936 865 790 085(43) 


2.977 274 273 728(33) 


4.984 713 354 070(45) 


6.989 892 949 082(32) 


1 


0.349 595 330 721(5i) 


2.758 891 177 876(36) 


4.851 946 642 761(42) 


6.900 301 395 128(3b) 


2 


-0.337 237 264 447(5i) 


2.487 025 791 777(3s) 


4.709 976 255 628(42) 


6.803 992 334 705(34) 


5 


-2.549 035 191 007(53) 


1.494 183 218 341(39) 


4.268 043 172 724(45) 


6.534 685 249 316(35) 


10 


-6.529 142 779 202(6o) 


-0.660 939 314 881(4o) 


3.318 493 978 272(46) 


6.100 400 048 017(38) 


12 


-8.182 546 155 166(6b) 


-1.659 292 230 771(44) 


2.838 014 627 229(48) 


5.905 881 549 211(3g) 


50 


-41.876 959 736 225(37) 


-26.863 072 307 493(33) 


-14.310 287 343 156(28) 


-4.206 192 073 796(3i) 



Acknowledgments 



Partial financial support of this work under Grant Nos. GP249507 and GP3438 from the Natural Sciences and 
Engineering Research Council of Canada is gratefully acknowledged by two of us (respectively NS and RLH). 



[1] J. F. Carincua, A. M. Perelomov, M. F. Raiiada and M. Santandcr, A tjuantum exactly solvable nonlinear oscillator related 

to the isotonic oscillator, J. Phys. A: Math. Theor. 41 (2008) 085301. 
[2] B. Champion, R. L. Hall, and N. Saad, Asymptotic Iteration Method for singular potentials. Int. J. Mod. Phys. A 23 

(2008) 1405 - 1415. 

[3] H. Ciftci, R. L. Hall, N. Saad, and E. Dogu, Physical applications of second-order linear differential equations that admit 

polynomial solutions, ,J. Phys. A: Math. Theor. 43 (2010) 415206. 
[4] H. Ciftci, R. L. Hall, and N. Saad, Asymptotic iteration method for eigenvalue problems, J. Phys. A: Math. Gen. 36 (2003) 

11807. 

[5] H. Ciftci, R. L. Hall, and N. Saad, Construction of exact solutions to eigenvalue problems by the asymptotic iteration 

method, J. Phys. A: Math. Gen. 38 (2005) 1147. 
[6] J. M. Fellows, R. A. Smith, Factorization solution of a family of quantum nonlinear oscillators, J. Phys. A: Math. Theor. 

42 (2009) 335303. 

[7] R. L. Hall, N. Saad, and A. B. von Kcviczky, Spiked harmonic oscillators, J. Math. Phys. 43 (2002) 94. 
[8] A. Ronveaux (ed.), Heun's Differential Equations, Oxford University Press, New York (1995). 
[9] J. Sesma, The generalized quantum isotonic oscillator, J. Phys. A: Math. Theor. 43 (2010) 185303. 
[10] N. M. Temme, Special functions: an introduction to the classical functions of mathematical physics, Wiley, New York, 
(1996). Laguerre polynomials are discussed in chapter 6. 



